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Abstract 

We discuss models for coupled wave equations describing interacting 
fields, focusing on the speed of travelling wave solutions. In particular, 
we propose a general mechanism for selecting and tuning the speed of 
the corresponding (multi-component) travelling wave solutions under cer¬ 
tain physical conditions. A number of physical models (molecular chains, 
coupled Josephson junctions, propagation of kinks in chains of adsorbed 
atoms and domain walls) are considered as examples. 


1 Introduction 

In a previous paper [1] dealing with a concrete physical model - more specifi¬ 
cally, with the nonlinear dynamics of the DNA macromolecule [2, 3] - we have 
observed a remarkable phenomenon. This is as follows: in the continuum limit 
that model reduces to two coupled nonlinear wave equations for different fields 
<^ 1>2 ; if the coupling is switched off, each of the wave equations £ 4,2 obeyed by 
the fields (j) 1,2 is Lorentz invariant with different limiting speed, i.e. in par¬ 
ticular admits travelling waves (solitons) with any speed c smaller than the 
limiting speed cy 2 . When the full model, including the interaction, is consid¬ 
ered, there is no Lorentz invariance, and the travelling wave solutions admits 
only a given speed (see also [4]). Thus we have a selection mechanism for the 
speed of travelling wave (TW) solutions (the latter turn out to be also stable 
and thus physically relevant [3, 4, 5, 6, 7]). 

In the present paper we want to study if this “speed selection mechanism” 
works also in a more general class of equations. We will answer to this in the 
positive, and actually the relevant class of equations turns out to be rather 
ample. An important point in this context is the fact that the speed selection 
mechanism of the paper mentioned above is implemented using a constraint [1], 
In general the constraint will define a submanifold of the configuration space 
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for the associated dynamical system. The crucial requirement for our speed 
selection mechanism to work is that the constraint is natural, i.e the associated 
submanifold in invariant under the dynamics. In other words, if the initial data 
are chosen to lie on this submanifold, the dynamics will take place entirely on 
it with no need to introduce external forces to enforce the constraint. 

We note that, apart from the theoretical interest, this question also have 
potentially very relevant practical consequences. In fact, for instance, it would 
point out a way to have transmission lines with physically determined speed for 
the travelling wave packets. 

As mentioned above, we will answer in the positive to the question of appli¬ 
cability of the speed selection mechanism to more general equations. It turns 
out that, albeit we are mostly interested in nonlinear equations, a speed se¬ 
lection mechanism is also present in the case of linear systems. Note that for 
linear equations the uncoupled equations have each a well definite speed, so in 
this context the mechanism we study amounts to a change in the value of the 
n allowed speeds. 

We will thus start, in Section 2, by considering simple linear systems; in 
this case the speed selection follows from some trivial algebraic facts and gives 
the possibility of tuning the speed of TW just by changing the value of an 
interaction parameter. 

The general setting of this note will be as follows. We investigate wave 
equations for N fields <f> 1 (x , t ), i = 1...N in 1 + 1 dimensions (one space di¬ 
mension and time) described by a Lagrangian, and we are specially interested 
in travelling wave solutions. 

The Lagrangian will be written as 



( 1 ) 


where the A is a Lagrangian for each of the fields </>*, and Ant is the interaction 
Lagrangian. 

We will make “minimal” choices for A and Ant, as we want to understand 
the phenomenon of speed selection in the simplest possible terms. In particular, 
Ant will be made of a “gradient interaction term”, coupling the spatial gradient 
(j) l x of different fields - and playing an essential role in our analysis - and possibly 
of a potential term V(</>*) = Vfy 1 ,..., 

We will consider mainly Lagrangians leading to hyperbolic wave equations. 
In appendix B, we will extend our considerations to parabolic, Schrodinger-like 
equations. Thus, we take first a Lagrangian of the form 


N N 



i=l ij =1 


where pi , tq are constant, 7 ^ are the components of a constant (N x N) matrix 
T, and the potential V (which could be zero) will be appropriately chosen below 
(see Sect. 2 and Sect. 3). The fields and the constants are assumed to be real; 


2 



note that the speed of the waves for the fields fa we get in the decoupled case 
7 = V = 0 are 

c» = Ki/pi . (3) 

The matrix T can be taken to be symmetric, and in the present notation 
we can assume it has zero terms on the diagonal, as the corresponding terms 
are represented by the Hi (so that T only represents the interaction between 
gradients of different fields). 

The Euler-Lagrange equations corresponding to the Lagrangian (2) are 1 

Pi <t>u - K? - 7 ij (fix = - {dV/dfa) . (4) 

It is appropriate to mention immediately some physical relevant cases in 
which one meets Lagrangians of the type (2). This include e.g., beside the 
“composite” model of DNA dynamics [2] mentioned above and the strictly re¬ 
lated case of long wavelength excitations in a chain of double pendulums (this 
also applies to polyethylene [ 8 ]), the case of coupled Josephson junctions [9, 10] 
and interaction between kinks in coupled chains of adsorbed atoms [11]. These 
cases will be discussed later on as examples of our general mechanism; see Sec¬ 
tions 4 and 5. 

Our analysis would of course also apply to Lagrangians differing from (2) 
only by boundary terms and total differentials, e.g by a term ^ JA ._ x aijfad x (fj', 
we focus on Lagrangians of the form (2) both for these make the analysis rather 
transparent and for their physical relevance. 

Most of our discussion will be conducted in the simplest nontrivial case, 
N = 2; we will then discuss the generalization to the arbitrary N case, which 
will be in some case rather immediate. 


2 The linear case 

As mentioned above, we start discussing the case of a quadratic Lagrangian, 
hence of linear field (Euler-Lagrange) equations. 

Moreover, we will at first consider the case TV = 2; we will then write, for 
ease of notation, fa = ip and fa = if', there is only one nontrivial term in the 
matrix T, i.e. 712 (which, with a slight abuse of notation, we denote simply by 
7 ). Thus we study the Lagrangian 

£ = \ [{Pi<Pt + pbli) - ( K i Pi + «l V>*)] - IPxfa - V(ip,ip) . (5) 

The Euler-Lagrange field equations are then 

Pi 'ptt 'pxx T "0 XX = {&V / dtp) 5 

plfat - kI^xx - IPXX = -{dV/dfa. (6) 

^^Here and below we move the field index up and down for typographical convenience. 
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In the simplest case, which we consider in this section, V will be a quadratic 
function of its arguments (including the case where V is trivial) and the equa¬ 
tions ( 6 ) will hence be linear. 

We will write, for the sake of definiteness, 

T(^,7) = \{p\v 2 + 2A<£■!/' + P 2 V’ 2 ) > (7) 


where p 2 2 are the masses of the non interacting fields (A = 7 = 0) and can 
hence be assumed to be positive. 

Thus our Lagrangian is defined by three matrices Q^\ Q^ x \ Q^ v \ given by 


Q (t) 



Q {x) 




g (y) 



( 8 ) 


we have (returning for a moment to the notation with indices) 


C = 


1 


N 

*J=1 


[(tiQflti) (^QY/o’r) (^gV) 


The problem - or better the source of interesting behavior - lies in that, 
in general, the three matrices Q ( ' t \ Q ^ and Q ^ do not commute with each 
other. More specifically, we have 


Q (t \Q (x) 

Q {x) ,Q {v) 

Q {t) ,Q {v) 


(pI pi) ( 7 0 7 ) ; 

(A-A) (" 7) + (4-A) (° 7) ; (9) 

(A-A) (° 7 ) ■ 


This lack of commutativity between the different matrices, and in particular 
the fact we have [Q^\Q^} 7 0 , is responsible for the breaking of the space-time 
Lorentz symmetry, which is one of the sources of the unusual (and interesting) 
features of the theory; see also the discussion in Appendix A. 


2.1 Travelling wave ansatz 

If we look for travelling wave solutions, i.e. for solutions of the form 

<p(x,t) = <p(x ± ct) = tp(z) , ip(x,t) = ip{x ± ct) = > (10) 

we have then to study the equations 

{pic 2 - nf) p zz -7 ip zz = - {p\<p + Mp) , 

{pic 2 - nf) ip zz - 7 ^ = - (X<p + /Z 2 VO ■ 
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( 11 ) 









These are immediately written in matrix form. Defining 


_ ([*l-Pic 2 ] 7 

' 7 [«1 -Pic 2 } 


B = 


rfD* = 


A 


we recast the previous equation as 

. 


(12) 


If det(A) 0, i.e. if 




zpipi 


(kIpI + 4pI) ±{k\pI-4p\) 


47 2 PiPi 


(«lP2 - «1 p?) 2 


we can invert A and further recast (12) as 


M 

M 

det(M) 


M $ , 

A- 1 B = M 


1 


/ 7 A - kIpI + c 2 p\pl 
\ — k 2 A + 7/t 2 + c 2 Apf 


7 2 -(^- C 2p 2 )(K 2 -c2p 2 ) 

Ad (A 2 — /i 2 //f) . 


—«|A + 7/z| + c 2 Apf \ _ 
7A - n\pl + c 2 pIp\ ) ’ 


One could then diagonalize the matrix M by a change of basis $ = A4> so 
that in the new coordinates the equation reads = M<f> with M = A -1 MA = 
diag(<5i(c), 6 2 {c)) a diagonal matrix, and now solve easily for $( 2 :) and hence for 
$( 2 ) = A$(,j). As the functions 8 i(c) (which of course do also depend on the 
various parameters of the equations) will in general satisfy <5i(c) y^ 82 (c), we will 
have two normal modes with frequencies (in z) satisfying (wj;~^) 2 = —<$i(c) 

and (u> 2 2 ^) 2 = — 82 (c). These relations correspond to the two branches of the 
dispersion relations we will find in the next section using the Fourier transform. 

Note that the situation is completely different in the V = 0 case, i.e. for 
B = 0. In this case the equation (12) reduces to 


A<f> zz = 0 , 


(13) 


which admits a solution if and only if det(A) = 0, which should be seen as a 
requirement on the speed c. More specifically, this yields with simple algebra. 


4 = 5 


1 r 


(c\ + c 2 2 ) ± (cl - c\) Vp 


(14) 


where we have written 


P = 1 + 4 


7 


_Plp 2 (c\ ~ 4 ). 

here and below c z = \ni/Pi\ is the speed of decoupled waves defined in Eq. (3). 
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The eigenvectors of the matrix A corresponding to the eigenvalues c± can 
be also easily computed; up to a normalization factor, they are 


$± = 


(ci /c 2 ) y/( ci-4)/( c i-4) 


(15) 


It is worth stressing that equations (14) and (15) describe two remarkable 
features of TWs in linear coupled systems (11) 

First, we have that the linear coupling between the <p and ijj field in the 
wave equation (11) forces a synchronization of the ip- and ip- waves: they have 
to propagate with the same speed, given either by c+ or by c_. To see this 
let us first diagonalize the kinetic part of the Lagrangian (5). It is not a priori 
evident that this is possible because the kinetic matrices Q x e Q l in Eq. (8) do 
not commute. However, one can show that a GL(2,M) transformation S exists, 
which acting on the vector $ as $ = S& diagonalizes the Lagrangian (5) with 

V = 0: 


C = 


1 

2 




(16) 


We see that in the $ frame we have two decoupled normal modes given by (15) 
propagating respectively at speed c+ and c_ given by Eq. (14); however they 
represent disjoint sectors, which cannot be superimposed. 

Second, we can tune the speeds of the tp and ijj waves between zero and 
a maximum value by changing the coupling between the two fields, which is 
parametrized by the (positive) ratio r, 

0 < r := (< c\ c\ := r max . (17) 

\P1P2J 

In particular, let us start from the decoupled case; here r = 0 and we get c+ = C\ 
and c_ = C 2 - If we now increase r, then the value of c_|_ increases whereas that of 
c_ decreases. When r —> r max , c+ tends to its maximum value whereas c_ —> 0. 
Notice that for r > r max , c 2 _ become negative, i.e we have an imaginary speed 
for the harmonic mode; this detects an instability of the system. Notice also 
that Eq. (17) defines a general bound on the coupling parameters, which can 
be also generalized to the case of a non vanishing potential (7). 


2.2 Fourier transform 

To solve the general linear case we well make use of the linear character of 
the field equations (6), and pass to consider the Fourier transforms f(q,ui) and 
g(q,oj) for, respectively, the fields <p(x,t) and xj)(x,t)', the equations (6) are then 
recast as 

- (p\uj 2 - K\q 2 ) / + 7 q 2 g + ( p 2 J + Xg) = 0, 

- {pW - i4q 2 ) g + 7 </ 2 7 + (A 7 +A< 2 ?) = ( 18 ) 
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These are promptly rewritten: defining now 


M 


([pW-kW-^1] - 7 g 2 -A \ p _ (f 

l -79 2 -A IpW-kW-hI}) ’ \g. 


equations (18) read 

MF = 0 . (19) 

Again we must require the vanishing of a determinant, i.e. det (M) = 0; this 
will now give a relation between uj 2 and q 2 , i.e. we will get some dispersion 
relations (DRs). 

More precisely, these read 


<4 = 


1 ■ 

2 . 


(c\ + cl)q 2 + (u\ + u\) ± VP 


( 20 ) 


where we have written 


P 


[/ 2 2\ 2 , / 2 2 \ 1 2 , a (79 2 + A) 2 

(ci-c 2 )g +{u 1 -u 2 )\ + 4- 2 2 — 

P 1 P 2 


with Ui = Hi I pi and c, as above. 

For a general nontrivial potential, we have two different DRs, which will also 
determine the phase velocity v± = dui±/dq of the wave. In this generic case the 
DR have a rather involved form. They take a simpler form for some particular 
or limiting cases. 

For a vanishing potential i.e pi = H 2 = A = 0, one can easily check that the 
DR (20) become linear (acoustic) and that duj±/dq = c± with c± given by (14). 
This case describes also the high energy limit q —> 00 of the DR (20), where the 
potential can be neglected. 

In the low-energy limit q —> 0 the DR (20) describe two optical branches, 
whose analytic expression can be derived expanding ( 20 ) near q = 0 . 

As expected, the DR become simple also in the decoupling limit 7 = A = 0. 
In this case, we have two branches with optical DR ui 2 2 = c 2 2 q 2 + u\ 2 . Another 
interesting particular case is c\ = C 2 , u\ = U 2 - We obtain also in this case two 
optical branches: 


<4 = [ci ± 


P1P2 


q 2 +u\± 


P 1 P 2 


( 21 ) 


The synchronization and tuning effect discussed in Sect. (2.1) for the TW 
speed c applies also to the phase speed du>±/dq. Now synchronization simply 
means that a given phase of the ip and ip wave must propagate at the same 
speed, whereas tuning means that we can change the phase velocity by acting 
on the coupling parameters of the two fields. In the case of acoustic and optical 
dispersion relations this tuning is very simple. In the acoustic case, the group 
and the phase speed are the same, and thus it has been already described in 
Sect. (2.1). In the optical case we have instead du>±/dq = c±q/uj ±, so that 
the phase velocity becomes zero whenever c± = 0 and for a given phase, grows 
monotonically with c±. 
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2.3 The N-fields case and Lorentz symmetry 

The discussion conducted above is readily generalized to the case of N fields. 
This is specially transparent in terms of the travelling wave ansatz (TWA): in 
this case we get again a matrix equation, but the matrix is now an (N x N) 
matrix. Solving the eigenvalue problem we get in general N normal modes. 

In terms of the Fourier transform, we will get N branches of the dispersion 
relations. 

Similarly to the N = 2 case, for a vanishing potential we will have N de¬ 
terminations of the allowed speeds c and we can change the values of the these 
speeds (for given m, pi and hence “uncoupled speeds” Cj = /?*/ pi) by acting on 
the coupling parameters 7 ^ of the model. 

Let us now briefly comment on the Lorentz symmetry of our two-dimensional 
field theory (2). Space-time Lorentz symmetry is explicitly broken in the La- 
grangian (2). This breaking has two sources. The first is the presence of several 
different limit speeds cp the breaking of the Lorentz symmetry in this case is 
expressed by the non commutativity of the matrices and QW in Eq. ( 8 ). 
The second source is the non covariance of the kinetic coupling term, which 
contains the space but not the time derivative. On the other hand, by imple¬ 
menting the TWA or by solving the field equations in Fourier space the relevant 
matrices are not the individual Q^ and but a linear combination of them. 

In Fourier space the breaking of the Lorentz symmetry is evident from the 
form of the dispersion relations (20). They are invariant under Lorentz boosts 
only when ci = C 2 and 7 = 0 . Despite of the explicit breaking of the Lorentz 
symmetry, it is quite evident that some remnant of it survives in the field theory 
(2). This is in particular evident in the massless case where the dispersion 
relation remain linear, the Lagrangian can be diagonalized describing decoupled 
sectors and the presence of the coupling term 7 just changes the values of the 
two speeds C 12 . We will discuss this remnant Lorentz symmetry and related 
group theoretical aspects in Appendix A. 


3 The nonlinear case 

Our discussion in Section 2 makes a substantial use of the linearity of field 
equations; thus it cannot be extended to the nonlinear case. 

In the nonlinear case, the TWA produces a system of nonlinear ODEs, and 
TW solutions are obtained as solutions cj>i(x±ct) of these nonlinear ODEs, which 
can be considered as describing an equivalent mechanical system; this procedure 
leaves (in general, see below) the TW speed completely undetermined. 

This is particularly evident when when we start from a Lagrangian which is 
Lorentz invariant: then the TW solution must be a function of the relativistic 
gamma factor, tfri = <fii((x±ct) / y/l — c/c) (where c is the limit speed). The wave 
speed c can be therefore arbitrarily changed in the range [0, c] by a Lorentz boost. 
Thus any speed selection mechanism in the nonlinear case must necessarily start 
from a Lagrangian in which the Lorenz invariance is explicitly broken to remove 



this degeneracy in c. 

Any speed selection mechanism must constrain the dynamics of the equiva¬ 
lent mechanical system to happen in a submanifold of the configuration space 
in which c is fixed. Moreover, if the submanifold corresponds to a natural con¬ 
straint , i.e. is invariant under the dynamics, the constraint can also be realized 
as a selection on the initial data. A simple realization of such a mechanism has 
been proposed in Ref. [1], Here we discuss in detail its dynamical implications, 
in particular existence of invariant submanifolds and the naturalness of the con¬ 
straints that can be used to select the speed of TW solutions. Our main goal is 
obviously the generalization of the results of Ref. [1]. 

We are now considering a Lagrangian of the form (2) with a generic (analytic) 
potential V(<j>i ), and Euler-Lagrange equations (4). In the case N = 2 these 
reduce to (5) and (6) respectively. 

Note that the kinetic part of the Lagrangian - and hence the second order 
terms in the Euler-Lagrange equations - is still characterized by the two matrices 
QW and Q^ x \ which in general do not commute. This non commutativity 
expresses the breaking of the space-time Lorentz symmetry at the Lagrangian 
level. On the other hand, as already noticed in the discussion of the linear case, 
after using the TWA (or passing to Fourier space) the relevant matrix will be 
a linear combination (with coefficients depending on the speed of the TW) of 
Q« and Q^ x \ 

With the travelling wave ansatz (10), we are led to consider the equations 

(pfc 2 -^)4, - 7 atiz = - VV/dP) . (22) 


3.1 The two-dimensional case 


We will again start by considering the case N = 2; with the notation introduced 
above for the linear case, we are thus dealing with the equations 


(«i - p\c 2 ) ip zz +7 i> Z z = ( dV/dip ) := /(^,- 0 ) , 
(«!- plc 2 )^ zz + 7 Vzz = (dV/dip) := . 

These are rewritten in matrix form as 


(23) 


M* zz = F , 


having of course written 


M = 

( ( K i - Pic 2 ) 

7 A 

. _ 

f mi 

1 ) 


V 7 

(«2 - p¥' ) J 


V 7 

m 2 J 



This matrix M combines the Q^ and QW 

Provided det(M) ^ 0, i.e. 7 2 ^ mim 2 , we can now consider a linear change 
of coordinates in the fields space, 


= rf 


9 


which diagonalizes (23). 2 

We can actually ask more, i.e. that A T MA = I. This is obtained e.g. for 

A = ( a _ -\/g(oQ/V m 2 {m 1 m 2 - y 2 )_\ 

\-[a"/+ ^/q(a)\/m 2 [ay 2 - mim 2 a ± 7 y / g(a)]/[m 2 ^/mim 2 - 7 2 / 

where a is a free parameter and 

q(a) := m 2 + (y 2 — mim 2 )a 2 . 

In terms of the fields rj , our equations (23) read now simply 

d 2 r]i dW 
~d. i 2 " = ~~ ’ 

thus the dynamics of the system is described by the motion of a particle of unit 
mass in the effective potential W(rf) := V[(f>(r])\. 

In dealing with PDEs one should also specify a function space for the search 
of solutions. In view of the physical meaning of the equations, it is natural to 
look for solutions of finite energy. In the case of a natural Lagrangian like (2), 
this means that the solutions ip(x , t) should go to a minimum of the potential 
for x —► ±oo; once we proceed to reduction to a system of ODEs via the TWA, 
this means that we look for solutions which go to a minimum of V (we stress 
the condition is on V, not on W) for z —> ±oo. 

It is easily seen that if minima of V correspond to minima of the effective 
potential W, these solutions are necessarily trivial. If, on the other hand, min¬ 
ima of V correspond to maxima of W, nontrivial solutions with the prescribed 
asymptotic behavior can exist; moreover they can either be doubly asymptotic 
to the same local minimum of V, or connect two distinct local minima, corre¬ 
sponding to the same energy level. 

Thus the existence of TW solutions with the relevant behavior depends on 
the sign of V/W, i.e. on the sign of the determinant of the Jacobian matrix 
A. This in turn will depend on the value of c 2 . (The discussion of concrete 
examples will better clarify this point.) 

Provided this condition is satisfied, we are thus searching for solutions to 
our equivalent mechanical system with prescribed limit conditions for r/(z ) and 
rj z (z ) at z —^ ±oo. In particular, we require 

lim 77 .( 2 ) = 0 ; lim r] z (z) = r]± 

z —>-±oo z—>± 00 

with rj± corresponding to (necessarily degenerate) local minima of the potential 

V. 

Such solutions are in general not unique (nor stable, even in the realm of 
solutions satisfying the same limit conditions), as immediately shown by the 

2 With V(rj) = V[<f>(r ))], the latter reads MAr/ zz = (A -1 ) T (dV/dr/), so that we have to 
look for A such that A T MA is diagonal. 
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example of a doubly periodic potential, i.e. of a dynamical system on the torus 
[12, 13, 14]. 

The situation - and hence our analysis - is greatly simplified if there are 
some one-dimensional invariant submanifolds connecting two degenerate min¬ 
ima. Needless to say, this is a non-generic feature, and in general the existence 
of such invariant submanifolds is related to the presence of some symmetry in 
the potential 3 . This could be a Noether symmetry, guaranteeing the presence 
of a conserved quantity and hence the reduction to a lower dimensional effective 
dynamics; or even a discrete symmetries - such as a reflection symmetry - in the 
potential. More generally it suffices to have a conserved or even a conditionally 
conserved quantity [15]; the existence of this is implied e.g. by the presence of 
a reflection symmetry. 

Albeit in general a conditionally conserved quantity and hence a low-dimen¬ 
sional invariant manifold is not necessarily associated to a symmetry, in practice 
this is most often the case. 

In this paper we will limit our considerations to the case in which the system 
has a reflection symmetry, so that we are guaranteed of the existence of a low¬ 
dimensional invariant submanifolds. We have now to distinguish between three 
possible cases: 

1) The system is fully invariant under Lorentz symmetry, i.e. we have c\ = 
C 2 = c and 7 = 0 in Eq. (23). In this case the invariant submanifold exists 
for every value of the soliton speed c below the limit speed c. We do not 
have selection of the soliton speed. The soliton speed is fully protected by 
the Lorentz symmetry. We will give examples of this situation in Sects. 
5.2 and 5.4. 

2) Lorentz symmetry is broken by 7 ^ 0 but we still have a single limit speed, 
i.e Ci = C 2 = c. Also in this case the invariant submanifold exists for every 
value of the soliton speed c below a limit speed and the soliton speed is 
protected by the Lorentz symmetry but the limit speed is changed by the 
kinetic interaction term. This is fully consistent with the discussion of the 
linear case of Sect. 2, where we have shown that the kinetic interaction 
term changes the wave speed according to Eq. (14). This situation is 
realized in an example discussed in Sect. 5.2. 

3) Lorentz symmetry is broken by ci C 2 - In this case the soliton speed is 
not protected by the Lorentz symmetry. This is the case in which we can 
have speed selection for the soliton; the invariant submanifold exists only 
for selected values of c. This situation is met in the examples discussed in 
Sections 5.1, 5.3 and 5.4. 

Summarizing, a sufficient condition for the speed selection mechanism con¬ 
sidered here to work is the presence of a reflection symmetry defining invariant 

3 Strictly speaking, we are here concerned with symmetries of the effective potential; these 
will however correspond to symmetries of the original physical potential, as the transformation 
mapping the latter into the former is smooth (and actually linear). 
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submanifolds connecting two degenerate minima of the potential W existing only 
for specific values of the soliton speed. Whereas a necessary condition for this 
speed selection mechanism to work is an explicit breaking of the Lorentz sym¬ 
metry at the Lagrangian level trough the presence of different limiting speeds. 


The practical implementation of the speed selection mechanism requires the 
determination of the invariant manifold associated with the reflection symme¬ 
try. This is a rather involved problem, which we will tackle using an ansatz. 
Obviously the manifold we are looking for must be invariant not only under the 
dynamics but also under the action of the reflection transformation. Because 
the kinetic part of the field equations (23) is linear the most natural ansatz for 
determining the invariant manifold is a linear equation involving the fields ip 
and ip. This will be our choice for the examples described in Sects. 5.1, 5.2, 5.3. 
It is also possible to search for invariant manifolds using non linear invariant 
functions involving the fields ip and ip. This will be our choice for the system 
described in Sect. 5.4. 

As a final step in the practical implementation of the speed selection mech¬ 
anism we need to restrict the dynamics on the invariant submanifold. This can 
be achieved in two different ways. A first way is through a constraint simply 
given by the ansatz we have used to determine the invariant manifold. Another 
way, making use of the fact we have a dynamically invariant manifold and hence 
a natural constraint , is to suitably select the initial conditions, i.e. choose them 
(position and velocity) along the chosen invariant manifold. 

3.2 The iV-dimensional case 

The analysis can be conducted along the same lines also in the general IV- 
dimensional case. Needless to say, even the preliminary step of diagonalizing 
the matrix M may be a substantial problem for high dimension; but we may in 
principles proceed along the same lines. 

Note that for dimension N a Noether symmetry will reduce the dynamics 
to a problem in dimension N — 1; thus we need N — 1 Noether symmetries in 
involution in order to reduce our problem to a one-dimensional one. Similarly, 
the presence of a reflection symmetry will in general only guarantee the existence 
of an invariant submanifold of dimension N — 1 (i.e. of codimension 1); thus we 
will need N — 1 reflection symmetries across planes with a nontrivial intersection 
to be sure of the existence of an invariant one-dimensional manifold. Note this 
is e.g. the case if the potential V{<pi) actually depends only on the squares 0? 
of the fields, or at least of N — 1 among them. 

Another situation which guarantees the existence of invariant one-dimensional 
submanifolds is that of a separable potential, V(0i,..., 0jv) = 14 (0i) + + 

Vn(<Pn), provided each 14(0fc)) has a minimum in 0*, = 0; in this case we do 
not have to require reflection symmetry of the 14 • 

One can of course also have a combination of the two situations mentioned 
above, i.e. a potential which is separable in potentials depending on several 
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groups of field variables, each of them having suitable (continuous or discrete) 
symmetries. 

Provided the potential has a sufficient degree of symmetry or separability, 
we are reduced to a one-dimensional analysis and we can proceed substantially 
as in the N = 2 case above. 

It is also possible that no symmetry is present, but that for specific values 
of the speed some one-dimensional invariant submanifold is present We expect 
however that this becomes increasingly unlikely - and anyway that it would 
be increasingly difficult to determine the allowed c in concrete terms - with 
increasing dimension N. 


4 Examples. Linear equations 

In this section we briefly apply our speed selection mechanism developed for the 
linear case to two examples. These describe the region of linear dynamics for 
some of the non linear systems to be considered in full in the next Section. 

4.1 Kinks in coupled chains of adsorbed atoms 

In the situation to be considered in Section 5.2, the linearized equations at 
(ip, ip) = (0, 0) read (setting /3 = 7 to comply with our nomenclature of coupling 
constants for the linear case) 

Ptt - Pxx = -tp - a (tp-ip) + 7 ipxx 

iptt - '‘Pxx = -ip - a {ip-ip) + 7 tp xx . (24) 

This linear system has exactly the form considered in Sect. 2, with special 
values of the parameters: 

Pi = P 2 = «i = K 2 = 1, pi = pi = 1 + a, A = —a (25) 

Here Ci = Ci = c = 1, Ui = thus Eq. (21) applies, giving the optical DRs 

= (1 ± 7 ) Q 2 + (1 + a) =F a , 

from which we get the phase speeds for the two branches 

C-(q) = — 7 = q ; c + (tf) = — . q . 

y/l + 2a+(l- 7 )q 2 v /1 + (! + 7 W 

Note that as a consequence of Eq. (17) we must have | 7 | < 1 and, as expected, 
for 7 = 1 (7 = — 1 ) only the w+ ( w_) branch is present, whereas the ui~ (w+) 
branch disappears. 

In the low-energy limit we get 

1 — ,_ y 

c -{q)---=== q , c+{q) ~ (1 + 7)9 ; 

Vl + 2a 

and in the high-energy limit 

c-(q) - V 1 ~ 7 , c+(q) ~ y/l + 7 • 
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4.2 Coupled Josephson junctions 

In the situation to be considered in Section 5.3, the linearized equations at 
(tp, ip) = ( 0 , 0 ) read (setting a = — 7 , again to comply with our general notation) 


Pit - Pxx = -<P + lip XX 

p 2 \pit - ipxX = - z' 2 ip + 7 ip xx . 


(26) 


This linear system has the form considered in Sect. 2 with special values of 
the parameters: 


Pi = «1 = «2 = 1, P‘2, = P, PI = 1, P2 = v, A = 0 (27) 

The dispersion relations (20) give 



The low-energy limit (q ~ 0) for the speeds c± = duj±/dq are given by 


c-(q) ~ — q + 0(q 3 ) ; c+{q) ~ q + 0(q 3 ) . 


In the high-energy limit (q 00 ) both speeds go to finite limits; the explicit 
expressions can be computed but are rather cumbersome and thus not reported. 

5 Examples. Nonlinear equations 

In this section we apply our TW speed selection mechanism developed for the 
non linear case to several examples. For all of these examples we look for an 
invariant submanifold of the dynamics for the associated mechanical system. 

5.1 Double pendulums chains 

The speed selection mechanism for TW solutions was originally proposed [1] in 
the context of mesoscopic models for non linear DNA torsion dynamics [2, 3, 4]. 
In particular, this was in studying the “composite Y model” [2] of DNA, in which 
the state of DNA is described by two angular variables ip(x, t) and d{x : t), which 
can be thought as describing (in the long wavelength limit, and thus using a 
continuum description) a chain of double pendulums. The peculiar geometry 
of DNA produces rather involved equations. The model is described by the 
Lagrangian density 


£ = ±{M 2 t -uH# 2 x + r 2 [mti 2 t + 


-u> s dl + 2(md t (tp t +d t ) - u s (ip x + d x )) cosip 

+m(ipt + ’dt ) 2 — u s (p x + 'dx) 2 ]} + 

+4r 2 A'p[cosi? + cos(y? + i?) — (1/2) COS 7 : — (3/2)] , (28) 
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and we refrain from writing the Euler-Lagrange equations for this Lagrangian; 
the reader is referred to [1] for the physical meaning of the various parameters 
appearing in C. 

Under the TWA, the Euler-Lagrange equations reduce to 

pp zz + /x(l+cos p)9 zz = 

= — 4K p sin(<p + 6) — /j sin(<p) ( 9 ') 2 + 2 K p sin(y>) ; 

/x (1 + cos p)p zz + [( J/r 2 ) + 2p(l + cosv?)] 0 ZZ = 

= -4K p (sin9 + sin(^ + 9)) + n sin(vj)[(^ 2 ) 2 + 2 p z 9 z ] , (29) 

where p := (me 2 — uj s ), J := ( Ic 2 — u>t), Wt = K t S 2 , ui s = K s 5 2 . 

The system in invariant under the reflection symmetry 9 — > —9, ip —> —<p. In 
the Lagrangian (28) appear two different limit speeds, thus space-time Lorentz 
symmetry is broken and we expect case 3) of Sect. 3.1 to apply, i.e that the 
invariant submanifold associated with the reflection symmetry exists only for 
fixed value of c. The natural candidates as invariant submanifolds of the dy¬ 
namics are therefore 9 = 0 and <p = 0. One can easily check that this is not the 
case for 9 = 0 , whereas setting ip = 0 in Eq. (29), and dividing the second by a 
factor 4, we get 

/i 9 ZZ = - 2 K p sin 9 , 0 ZZ = -2K p sm9 . (30) 

Thus ip = 0 is an invariant manifold if and only if J = 0; this fixes the speed of 
the TW to c = u t /1. 


5.2 Kinks in coupled chains of adsorbed atoms 

Let us now consider chains of adsorbed atoms (also called “adatoms”), and in 
particular the interaction between kinks in such chains [11, 16]. In our notation, 
these are described by the equations [11] 

<Ptt - Vxx = - sin ip - a sin(v? - ip) + /? ip xx , 

iptt - ipxx = -sin ip - a sin (ip - p) + Pp xx . (31) 


(These are also studied in [17] in the case /? = 0.) These correspond to the 
Lagrangian 


L = \ [(<fit - <P X ) + (V>? ^ V> 2 )] - P <Px tpx + [cos(^) + cos (VO + cos(<^ 

The TWA produces the equations 

(c 2 - 1 )ip Z z = ftipzz - sin p - a sin(y> - ip) 

(c 2 -l)ip zz = f3p zz - sin ip - a sin (ip - tp) , 


VO] • 

(32) 


which are written in the form 


M $" = F 
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by setting 


M = fM • F = MU “ asinl ! v ~' l, \) ; 

Y —p c — 1 y \ — smip — asm (ip — ip) J 

note that det(M) = (c 2 — l) 2 — P 2 ■ Under the condition c 2 ^ l±/3, the equation 
is rewritten as 

(1 — c 2 ) sin</3 + a(l + ft — c 2 ) sin(ip — ijj) — psimp 
= /3 2 - (1-c 2 ) 2 ’ 

( (1 — c 2 ) siny> + a(l + P — c 2 ) sin(ip — tp) — (3simp 

* = P 2 - (1-c 2 ) 2 ' 

This obviously admits two discrete symmetries: 

(tp,ip) ^ (ip,<p) and (p,ip) ->■ (-ip, -ip) . 


The latter is of no use (the invariant set it identifies is just ip = ip = 0, which 
gives the trivial solution), while the former suggest to pass to field coordinates 


p - ip p + ip 

V ■= —g— , £ “ — 

In terms of these, the equations are rewritten as 


V ZZ = 


asin(2ry) + sin(ry) cos(£) 2 a cos(77) + cos(£) 


£ zz = - sin(£) 


(c 2 - 1 )+P 
cos(r]) 

(c 2 - 1 )-p ■ 


(c 2 -l) + , 


We see immediately that both the submanifolds identified by £ = 0 and by 77 = 0 
are invariant. 

In the equation of motion (31) appears one single limit speeds, which with 
the units used in (31) is Ci = C 2 = c = 1, and a non vanishing kinetic coupling 
term p, with |/3| < 1. Thus space-time Lorentz symmetry is broken only by this 
coupling term, and case 2) of Sect. 3.1 applies. We therefore expect that the 
invariant submanifold associated with the discrete symmetry exists for every 
value of c below a limit value depending on p. We will now show that this is 
indeed the case. 

For P = 0 we would have a motion in an effective potential 


W = 


^2-- [cos(£ + 77) + cos(£ - 77) + a cos(2?7)] ; 


the analysis would be rather simple. The case P = 0 corresponds to not having 
the kinetic interaction term, see (31), (32). 

By restricting to 77 = 0 and writing hp. = [c 2 — 1 — /3 ] _1 we obtain 


P.ZZ = - h+ sin(£) , 


(33) 
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i.e. a standard sine-Gordon equation, which will support the standard sine- 
Gordon solitons [18, 19]. 

Note however that we have < 0 if and only if c 2 < 1 + /3; thus £- 
waves are possible if and only if the speed |c| is smaller than a critical value, 
|c| < Cf = y/1 + f3. 

By restricting instead to £ = 0 and writing now = [(c 2 — 1) +/3]” 1 we get 


r]zz = —h- sin( 77 ) [1 + 2a cos(? 7 )] 


(34) 


i.e. a double sine-Gordon equation, which will also support solitons [20]. 

Note that here too h- < 0 is not automatic: this is the case if and only if 
c 2 < 1 — (3, i.e. again 77 -waves are possible if and only if the speed is smaller 
than a limit value, which is now |c| < c v = y/1 — 0. 

Consistently with our discussion of the linear case, in order to have a real 
speed c we must require |/3| < 1 (corresponding to the condition |'y| < 1 used in 
Sect. 4). 

It should also be noted that for f} = 0 space-time Lorentz invariance if fully 
preserved, i.e it is realized the case 1) described in Sect. 3.1. The invariant 
submanifold exists for every value of the soliton speed smaller than c = 1 . 

5.3 Coupled Josephson junctions 

A weakly coupled system of two long Josephson junctions is described, in suit¬ 
able units, by the equations [9] 


<fitt ~ <Pxx = -sin(^) - a ip xx - erupt + fi , 

fJ 2 4>tt - ipxx = -v 2 sin(ip) - a tp xx - a 2 fJ.vip t + v 2 h ! 


we refer to [9] for the physical meaning of the different constants. Note that 
here cr, describe dissipation effects, and fi are external forcing needed to coun¬ 
terbalance dissipation. Considering the idealized case where dissipation is zero 
(and hence setting to zero also the external forcing terms fi), we are reduced to 


Vtt - Vxx = — sin(<£>) - a , 

1 i 2 iptt - ipxx = -v 2 sin(ip) - aip xx 


(35) 


This corresponds to a potential V = — [cos(v?) + v 1 cos(V’)]. 
The TWA reduces (35) to 


(c 2 -l)ip zz + a ip zz = -sin(v?) , 
(/r 2 c 2 — 1) ip zz + aip zz = -v 2 sm(ip) 


( 36 ) 


This is written as M<& zz = F if we define 
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Provided det(M) := [(c 2 — l)(/x 2 c 2 — 1) — a 2 ] ^ 0, and writing then h = 
[det(M)] -1 , the equations (36) are rewritten as <f> 22 = M~ 1 F, i.e. as 

<Pzz = — h [(1 - p 2 c 2 ) sin(y>) + a v 2 sin(^)] 

4>zz = —h [v 2 (1 - c 2 ) sin(ip) + a sin(t^)] . (37) 

We would like to determine one-dimensional invariant submanifolds for the 
dynamics (37). We look for discrete symmetries of the model. The system 
(35) is invariant under the reflection symmetry <p —> —ip, ip —> —ip. The most 
natural candidates for invariant submanifolds are therefore ip = 0 and ip = 0. 
These are, however, of no use, because correspond to trivial solutions of the field 
equations. Thus in general the determination of the invariant submanifold is a 
hard problem, but we can look for linear invariant submanifolds. 

That is we will look for linear combinations of the fields, £ = di<p + d 2 ip and 
?7 = d 3 tp + diip, such that the submanifolds £ = 0 and rj = 0 are invariant under 
(37). The inverse transformation can be written as 4 

<p = fci rj + k 2 i , ip = k 3 r] + fc 4 f ; 


in this way, and setting 

H = , , H , , = h p , 

AC1/C3 — /C2&4 

the equations (37) are written as 

izz = -H [[(afci — fc 4 (l — c 2 p 2 )] sin(fcir; + k 2 (p) 

+ [(1 — c 2 )fci — ak^v 2 sin(fc 4 ?7 + fc 3 £)] 

Vzz = -H [[(l-c 2 pi 2 )k 3 -ak 2 ]sm(k 1 r] + k 2 ^) 

- [(1 - c 2 )k 2 + ak 3 \v 2 sin(fc 4 77 + fc 3 £)] . 

On the line 77 = 0 these reduce to 

P.zz = -H [[(afci — fc 4 (l — c 2 /i 2 )] sin(fc 2 f) + [(1 - c 2 )fci - ak^u 2 sin(fc 3 £)] 
Vzz = -H [[(1 - c 2 p 2 )k 3 - ak 2 ] sin(fc 2 £) - [(1 - c 2 )k 2 + ak 3 ]v 2 sin(fc 3 £)] ; 

thus the line 77 = 0 is invariant if and only if 

(1 — c 2 p 2 )k 3 — ak 2 = 0 and (1 — c 2 )k 2 + ak 3 = 0. (38) 


This can be cast in the form of a matrix equation = 0, where we have 

defined 




l(l-c 2 ) a ) 


Ki 



4 Writing p = (k\k 3 — k 2 k^) 1 , the relation between the d; and the k, coefficients is given 
by d\ = —fc 4 p, d 2 = kip, d 3 = k 3 p, d 4 = —k 2 p. 
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The solution exists provided det(Q^) = 0, and this condition considered as an 
equation for c provides 


c' 


,2 



(l + /x 2 ) ± i/(l — / 2 2 ) 2 + Aa 2 p 2 

2~j? 


Notice that, in accordance with our results in the linear case, in order to 
have a real speed c we must require |a| < 1 (corresponding to the condition 
| 7 | < 1 used in Sect. 4). 

Denoting by Q({,±) the matrix Q 5 with c 2 = c 2 ^ ± ,, the kernel of Q(^,±) is 
spanned respectively by 


(1 — p 2 ) ± \/(l — p 2 ) 2 + 4 a 2 p 2 




2a 


One can check that with these choices for c and for Aq, Aq the line rj = 0 is 
invariant, and the evolution of £ is governed by an equation of the form 


t,zz = Pi sin(2o:£) + p 2 sin(u; ? , ± £) , 


where pi are coefficients depending on (hi, k 4 ), whose explicit expression can be 
easily computed but is not interesting here, and 


W( € .±) = -(1 + V) ± Vi 1 - M 2 ) 2 +4«V . 


One could analyze in the same way the invariance of the line £ = 0. Note 
however that there is no physical difference between the two fields, as they are 
just generic linear combinations of <p and ip-, thus the £ = 0 case will reproduce 
- with an exchange of roles between £ and 77 - the same results obtained for 
rj = 0. 

Note that in the case analyzed here the Lorentz symmetry is broken by the 
presence of two different limit speeds in the field equations (35). Thus, case 
3) of Sect. 3.1 applies. The invariant submanifold associated with the discrete 
symmetry exists only for selected values of the speed c. 

5.4 Modified Katsura model 

Our last example is a “two speeds of light'’ generalization of the model of Kat¬ 
sura [21]; this has been proposed to describe the coupling of magnetic and 
ferroelectric domain walls [21, 22]. The model we consider is described by the 
Lagrangian (for simplicity we set in the model of Ref. [21] the coupling constant 


A = H = 1) 



The Lagrangian has the discrete symmetry ip —> —ip, so that the natural 
candidate for the invariant submanifold for the dynamical system one obtains 
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after the TWA is ip = 0. Indeed such invariant manifold exists for every value of 
the parameters and the reduced dynamics is described by ip zz = ip 3 — 1. However, 
this equation does not support solitonic solutions because the potential does not 
have the required minima structure. 

The model also has another reflection symmetry, described by 

~ ■ 

If we look for analytic functions invariant under such a reflection, they are 
necessarily built as (algebraic) functions of the basic invariants 

ip 2 , cos 2 (/3^) , p cos/3^ . 

In fact, another invariant manifold for the dynamics can be found using the 
results of Ref. [21], where exact kink solutions have been found for /3 2 = 1/2 
using the ansatz 

ip = cos Pip. (40) 

One can easily show that Eq. (40) with f3 2 ^ 1/2 determines an invariant 
manifold for the equivalent mechanical system describing the TW dynamics 
stemming from (39), if the TW speed is fixed by 


c 


C2C1 


I 1-2/3 2 

4 - 2 / 3*4 


(41) 


Note that again we have two limiting speeds in the Lagrangian (39); thus the 
model falls in case 3) of section 3.1, and we do indeed have the behavior predicted 
there. On the other hand, in the case C\ = c 2 and f3 2 = 1/2, consistently with 
the results of Katsura [21], equation (40) describe an invariant manifold of the 
dynamics for every value of c; thus in this special case the system falls in case 
1) of the classification given in Sect. 3.1. 


6 Discussion. Field theoretical considerations 

This paper focused on the description of TW dynamics in macroscopic and 
mesoscopic systems. This is particular evident in the examples we have chosen 
in Sects. 4 and 5 to elucidate our mechanism. They represent macroscopic or 
mesoscopic models for molecular, biological or condensed matter systems. 

Nevertheless, at least some of the results of these paper have a much broader 
relevance and can be applied, generically, to any two-dimensional (2D) field 
theory. This is particularly true for the first part of paper, which concerned 
linear dynamics. 

An important point in this context is the presence in the action for the 2D 
field theory (2) of N different “speeds of light” Cj, see Eq. (3). This is a nat¬ 
ural condition for macroscopic and mesoscopic systems, where different modes 
propagating in a medium experience different effective physical parameters pi , Hi 
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(elastic, optical, magnetic and so on) of the media so that at the linear level their 
perturbations propagate at different speeds. On the other hand a fundamental, 
microscopic, 2D held theory has to be Lorentz invariant, thus characterized by 
a single speed of light. 

The case of a single limit speed is a particular case of the description of 
linear waves given in the first part of this paper (see e.g Eq. (21)). Thus, our 
results including those related to synchronization and tuning of the wave speed 
hold true for c\ = C 2 . 

From the Held theoretical point of view a particularly important case is 
that of a non interacting theory. When V = 7 = 0 and Cj = c the theory 
(2) is a 2D conformal Held theory (CFT) . It describes CFTs with N scalar 
fields (i.e. CFTs with central charge N), which play a fundamental role in 
several contexts such as string theory [23], critical points of phase transitions 
[24], microscopic explanation of black hole entropy [25], just to mention only 
few of them. The conformal invariance of the Lagrangian (2) is fully evident by 
passing to light cone coordinates x± = x ± ct. In this coordinates the theory 
is explicitly invariant under the action of 2D-diffeomorphisms (the conformal 
group in 2D) and the field equations read d x+ d x _ip = 0, whose general solution 
is a generic combination of right and left moving TW <p = /( x+) + f(x-). The 
conformal invariance is a general feature of the massless case. In fact it is also 
evident from Eq. (13), which does not fix the dependence on z of the fields 
ip(z), tp(z) but just gives a linear relation between the two fields. 

As a final remark we note an analogy between the synchronization mecha¬ 
nism for TW in the linear regime we have found in the first part of this paper, 
and Quantum Mechanics; this goes as follows. 

Let us consider for simplicity the field equations (11) with vanishing potential 
(our considerations can be easily extended to the V ^ 0 case). The system can 
be described in two different frames of the field space ip, 1 ft: a frame in which the 
two fields decouple completely i.e. in which the kinetic matrices are diagonal; 
and a frame in which the two fields interact with the interaction term 7 . In 
the first frame TW for the two fields can propagate independently, hence with 
different speeds c+ and c_ given by Eq. (14). Owing to the interaction, in 
the second frame the tp— and if)- waves are forced to travel with the same speed 
(synchronization). 

This relationship between TW waves in the two frames bears a strong anal¬ 
ogy with the eigenstates of two non-commuting quantum mechanical operators, 
acting on a 2D Hilbert spaces and associated to two non compatible physical 
observables. 


7 Conclusions 

In previous work, dealing with a concrete model for the nonlinear dynamics of 
DNA [2], we observed a peculiar mechanism which fixed the speed of solitons [1], 
In this paper we investigated if this mechanism could extend to a more general 
class of equations, answering this question in the positive. More generally, we 
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have considered systems of coupled wave equations; when decoupled each of 
them is Lorentz invariant with a limit speed Cj, with possibly different limit 
speeds and we have studied how the speed of travelling waves is affected by the 
coupling. 

In particular we have considered systems admitting a variational (Lagrangian) 
description, and in which the coupling between the different equations could in¬ 
volve a kinetic term, see Eq. (2). 

We observed that the general Lagrangian under consideration here, (see 
Eq. (2)) could fall in different classes according to its properties under Lorentz 
transformations. That is, Lorentz invariance could be unbroken (case 1), in 
which case there is of course no speed selection mechanism; it could be broken, 
albeit the coupled equations admit the same limit speed, by coupling terms 
(case 2), in which case there is still no speed selection but the limit speed can 
be changed due to the coupling; or finally the Lorentz invariance can be broken 
by the presence of different limit speeds and kinetic coupling terms (case 3), 
in which case we have a full selection of the speed of travelling wave solutions, 
which can only take a finite - and rather small, being limited by the number of 
coupled equations - set of values. 

In the latter - and more interesting - case, simple travelling wave solutions 
can be described as dynamically invariant one-dimensional manifolds for a me¬ 
chanical system (parametrically dependent on the speed c of travelling waves) 
associated to the system of PDEs under study, joining two extremal points for 
an effective potential which comply with limit conditions dictated by the finite 
energy condition for the PDEs system. Such manifolds only exist for special 
values of c, and this ignites the speed selection mechanism. 

We have then validated our general discussion by a number of physically 
significant examples; in particular, in Sects. 5.1, 5.3 and 5.4 we have shown that 
the speed selection mechanism here considered is present in double pendulum 
chains, coupled Josephson junctions, and in a modified Katsura model for the 
coupling of magnetic and ferroelectric domain walls. 

It is interesting that the mechanism described here - and which has a partial 
counterpart for linear equations as well - has some points of contact with general 
field-theoretic questions, as discussed in Sect. 6. 

In the Appendices, we give a closer look at some group-theoretical questions 
(Appendix A); and argue that, albeit here we worked specifically with hyperbolic 
PDEs, the same mechanism can be at work in some type of parabolic equations 
(Appendix B). 
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A Lorentz symmetry and group theoretical con¬ 
siderations 


In Sect. 2.3 we have noted that in the Lagrangian (2) the usual space-time 
Lorentz symmetry is explicitly broken by the presence of several limiting speeds 
Ci. However, there is some remnant of the Lorentz invariance; this can be 
described as follows. 

The Lorentz transformation 


x —> x 


{x — vt) 

y/1 -V 2 /C 2 ’ 


t —y t 


(t — ( v/c 2 )x ) 

y/l-V 2 /# 


is generated by the vector held 

X = 


x d 
c 2 dt 



(42) 


The evolutionary representative [26, 27, 28, 29] for this is the (generalized) vector 
held 

„ (x \ d 

Xev = y ^ ~b ^ ^ x J , 

where u is the dependent variable (the physical held). 

Thus if we consider the Lorentz symmetry with several limit speeds Ci acting 
on the held (j/, this is generated in the evolutionary representation by 


XI 



d 

d<p a 


If now we consider a generalized vector held 


X L 


XI 


X. 


N 


which (unless cf = ... = c 2 N ) will not be the evolutionary representative of any 
Lie-point vector held acting on the space of the (x, t; (j) 1 , ..., <j> N ) variables, this 
will be a symmetry for the Lagrangian C\ +... + Cn, where we dehne the partial 
Lagrangians as 

A = \ [Pi W4) 2 - «?(4) 2 ] . 

and for the corresponding held equations (j) l tt — c 2 (f> l xx = 0 as well . 5 

Note however that this will not leave the interaction Lagrangian, nor the 
corresponding interaction terms in the held equations, invariant. 

The symmetry properties of the Lagrangian and held equations described 
above can be also understood in terms of the commutation relations (9) of the 
matrices Q^\ Q^ x ' dehned in ( 8 ). We have a symmetry of the Lagrangian when 

5 Actually, here only derivatives of first (for the Lagrangian) and second (for the field equa¬ 
tions) order matter; so we could as well consider only transformations of the field derivatives, 
e.g. prescribing these are undergoing a hyperbolic rotation irrespective of any transformation 
on the field themselves. In this way one would be led to consider “hidden symmetries” [30]. 
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the matrices can be diagonalized simultaneously, i.e when 7 = 0 (corresponding 
to absence of coupling in the kinetic sector) and c* = c, i.e the same limit for 
the TW of different fields. 

The TV-held components TW will belong to the representation of the Lorentz 
symmetry generated by the generalized vector field Xl. More specifically they 
will transform via reducible representations 

T = Ti © T 2 © ... © T n , 

where each T a is a two-dimensional representation made of the Lorentz group. 
The TW solutions of the model with vanishing, respectively non-vanishing, po¬ 
tential belong to massless, respectively massive, representations of the Lorenz 
group. 


B Extension to parabolic equations 

We have so far considered coupled hyperbolic equations; on the other hand 
the treatment of the linear case is based on basic algebraic facts and does not 
involve hyperbolicity; it thus appears that our discussion can be extended to 
more general (time-evolution) equations, and in particular to parabolic ones. 

Here, we will not discuss the general case but only consider, for the sake of 
concreteness, TW dynamics described by coupled linear Schrodinger equations 
(LSE); this can be thought as emerging from the linearization of the nonlinear 
Schrodinger equations (NLS) 

iipt = + reM 2 '*/’ • (43) 

This is known to describe, among others, optical solitons and, upon introducing 
a potential term V(x)ip, Bose-Einstein condensates. As usual we will consider 
the N = 2 case, the generalization to the IV-helds case being performed along 
the lines described in Sect. 2.3. 

Let us consider the system of two coupled LSE, 


i 

1 

7 , 

dV 


— <Pt + 

Cl 

7 ^ ( !Pxx 

~2 rXX — 

~~d^> ’ 


i , 

1 

7 

dV 

(44) 

—i>t + 
C2 

2 rxx 

~2 ^XX 

~d^ ’ 


where c», 7 are some parameters and now the fields tp, tp are generically complex. 
In analogy with the hyperbolic case we have introduced a quadratic potential 
of the form given by Eq. (7). As in the hyperbolic case, to solve this equations 
we pass to the Fourier transforms. 

In the decoupled case 7 = A = 0 we get two dispersion relations 07 (g), i = 
1 , 2 and the related phase velocities v= du±/dq, 

Ui = ^Ci (g 2 - /z 2 ) , v {p) = dq ; (45) 
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thus Ci characterizes the phase velocity of the wave packet. 

In the generic coupled case 7 ^ 0, V ^ 0, we get dispersion relations similar 
to ( 20 ), i.e. 


W±= 2 


-(ci + c 2 )q 2 - (ci/i? + C 2 H 2 ) ± V'PiQ 2 ) 


(46) 


V(q 2 ) = ^(ci + c 2 ) V - (ci/i? + c 2 /r|)(ci + c 2 )g 2 + 

+(ciMi + ^2M2 ) 2 - cic 2 (<7 2 - 2 nl)(q 2 - 2 fj%) + CiC 2 (<7 2 7 + A) 2 . ( 47 ) 


Notice also that in the case of vanishing potential we get the dispersion 
relations 

w± = i c±q 2 , (48) 

where c± has exactly the same form (14) found in the hyperbolic case with 
squared-velocities replaced by the parameters Cj : 


c± 


1 

2 


Cl + c 2 ± (ci 


C2) \ /l + 


4 7 2 Ci c 2 
(ci - c 2 ) 2 


(49) 


As in the hyperbolic case this equation tells us that we can tune the phase 
velocity of the ip and ip waves just by changing the parameter |q| in the range 
[0,1], (see Eq. (17)). In particular, the phase velocity of a given phase in w+ 
(u>-) changes from a maximum (minimum) value when 7 = 0 , given by the 
smallest (greatest) of c±, to zero (a maximum value) when |y| = 1 . 
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